Flux-dependent occupations and occupation difference in geometrically symmetric 
and energy degenerate double-dot Aharonov-Bohm interferometers 
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We study the steady-state characteristics and the transient behavior of the nonequilibrium double- 
dot Aharonov-Bohm interferometer using analytical tools and numerical simulations. Our simple 
setup includes noninteracting degenerate quantum dots that are coupled to two biased metallic leads 
at the same strength. A magnetic flux $ is piercing the setup perpendicularly. As we tune the degen- 
erate dots energies away from the symmetric point we observe four nontrivial magnetic flux control 
effects: (i) flux dependency of the dots occupation, (ii) magnetic flux induced occupation difference 
between the dots, at degeneracy, (iii) the effect of "phase-localization" of the dots coherence holds 
only at the symmetric point, while in general both real and imaginary parts of the coherence are 
nonzero, and (iv) coherent evolution survives even when the dephasing strength, introduced into 
our model using Biittiker probe, is large and comparable to the dots energies and the bias voltage. 
Moreover, not only finite dephasing strength does not destroy the coherence features, it can provide 
new type of coherent oscillations. These four phenomena take place when the dots energies are 
gated, to be positioned away from the symmetric point, demonstrating that the combination of bias 
voltage, magnetic flux and gating field, can provide delicate controllability over the occupation of 
each of the quantum dots, and their coherence. 

PACS numbers: 73.23.-b,85.65.+h,73.63.-b 
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I. INTRODUCTION 

The steady-state properties of the Aharonov- 
Bohm (AB) interferometer have been intensively 
investigated^, with the motivation to explore coherence 
effects in electron transmission within mesoscopic and 
nanoscale structures^. Particularly, the role of electron- 
electron (e-e) interaction effects in AB intcrferomctry 
has been considered in Refs£~— , revealing, e.g., asym- 
metric interference patterns^ and the enhancement^ 
or elimination^ 2 - of the Kondo physics. Recent works 
further considered the possibility of magnetic-field 
control in molecular transport junctions^—. The 
real-time dynamics of AB interferometers has been of 
recent interest, motivated by the challenge to under- 
stand quantum dynamics, particularly decoherence and 
dissipation, in open nonequilibrium quantum systems. 
Studies of electron dynamics in double-dot AB inter- 
ferometers in the absence of e-e interactions have been 
carried out in RefsJ^— , using a non-markovian master 
equation approach. The role of e-e repulsion effects 
on the dots dynamics was studied numerically using a 
non-perturbative method in Refi^S. 

In this paper, we focus on a simple-minimal model, 
the noninteracting double quantum dot AB interferom- 
eter, and study its transient and steady-state properties 
in biased situations. For a scheme of this model, see Fig. 
[TJ This system has revealed a wealth of intricate be- 
havior, such as "flux-dependent level attraction"—, and 
the ability to achieve decoherence control when junction 
asymmetry is incorporated^. Here, this noninteracting 
system further displays other non-trivial effects in both 





FIG. 1: Scheme of a double-dot AB interferometer. The two 
dots are each represented by a single electronic level, which do 
not directly couple. The total magnetic flux is denoted by $. 
Dot '1' may be susceptible to dephasing effects, introduced 
here through the coupling of this dot to a dephasing probe, 
the terminal P. The role of dephasing effects is discussed in 
Sec. V. 



the transient regime and the stationary limit. While pre- 
vious studies have allowed for junction asymmetry and 
nondegenerate dot o 18 ' 21 , we restrict ourselves to the sim- 
plest case of energy degenerate dots and symmetric dots- 
lead couplings. However, in our study the degenerate 
levels may be tuned away from the symmetric point, a 
situation that can be reached by applying a gating field. 
Using exact analytic expressions and numerical simula- 
tions, we expose several nontrivial effects emerging in this 
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biased-gated AB setup: (i) First, dots' occupations dis- 
play strong flux dependency, (ii) Second, not only do the 
occupations vary with flux, but the dots acquire unequal 
occupations, at degeneracy, (iii) Further, we show that 
the effect of "phase localization" 17 appears only at the 
symmetric point, while when the system is gated away 
from that point dot coherences can be feasibly controlled 
by the bias voltage. 

In the stationary limit we further study the role of 
dephasing effects, implemented here through a dephas- 
ing probe^r— , on the coherence properties of the sys- 
tem. Interestingly, we find that finite and substantial 
dephasing strength, at the order of the bias voltage and 
dots energies, still allows for flux dependency of occupa- 
tion. However, at finite dephasing the occupation-phase 
dependency significantly differs from the zero dephasing 
limit. In other words, finite dephasing alter coherent os- 
cillations, to provide new features. The flux dependency 
of the occupation is fully washed out at a very strong 
dephasing strength. 

The structure of the paper is as follows. In Sec. II we 
present the double-dot AB interferometer model. Sec. 
Ill explores the steady-state properties of the system us- 
ing the nonequilibrium Green's function approach. We 
derive closed analytic expressions (at zero temperature) 
for the dots' level occupation, the coherence between the 
dots, and the charge current in the system. Sec. IV 
provides numerical results in the transient regime, indi- 
cating on the time scale it takes for the system to reach 
the stationary limit, and on the intricate dynamics in- 
volved. Sec. V. details the role of dephasing effects, 
providing analytic expressions for the dots occupations 
in the steady-state limit. Sec. VI summarizes our main 
results and concludes. 



II. MODEL 

We focus on a symmetric AB setup, with a quantum 
dot located at each arm of the interferometer. The dots 
are connected to two metal leads (refereed to as baths, or 
reservoirs) maintained in a biased state. For simplicity, 
we neglect the spin degree of freedom and describe each 
quantum dot by a spinless electronic level, see Fig. Q] 
The total Hamiltonian H includes the following terms 

H = H S + H B +H SB . (1) 

Here H$ is a subsystem (dots) Hamiltonian, H B includes 
the two metals and Hsb incorporates subsystem-bath in- 
teraction terms. Specifically, we assume uncoupled dots, 

H s — Eia[ai + e 2 a} 2 a 2 . (2) 

To keep our discussion general, we allow the states to be 
nondegenerate at this point. In our analytic and numeri- 
cal calculations below we have forced degeneracy, at and 
ap are the subsystem creation and annihilation operators, 



respectively, where (3 = 1, 2 denotes dots T' and '2'. The 
metals are composed of noninteracting electrons, 

H B = } J U3ia\ai + } j u r ata r , (3) 

l r 

where a\ r and a^ r are bath creation and annihilation 
operators, for the left (I £ L) and right (r g R) leads. 
The subsystem-bath interaction term is given by 

H SB = ^ a \ a ^ + E Cp,T4«f> ei *> + h - c -> ( 4 ) 
p,l p,r 

where £ is the coupling strength to left bath and similarly 
C stands for the coupling strength to the right bath. The 
notation here is general, but we later take these couplings 
to be identical, since we are interested in a dot-lead sym- 
metric setup. Here 4>^ and <j>^ are the AB phase factors, 
acquired by electron waves in a magnetic field perpendic- 
ular to the device plane. These phases are constrained 
to satisfy the following relation 

4>1 - 02 + <f>? ~ 4>2 = <f> = 27T$/$ . (5) 

$ is the magnetic flux enclosed by the ring and $o = he/ e 
is the flux quantum. In what follows we adopt the gauge 
<ft - 4% = <f>f- = 0/2. 

We voltage-bias the system, using the convention 
A/i = fiL — fiR > 0, with [j.l,r as the chemical po- 
tential of the metals. While we bias the system in a 
symmetric manner, /i£ = —(J,r, the dots levels may be 
placed away from the so called "symmetric point" at 
which (iL — £p = £p~ M-R- This situation may be achieved 
by applying a gate voltage to each dot. For simplicity, we 
use the following conventions below, h = 1 and electron 
charge e = 1. 



III. STATIONARY BEHAVIOR 

A. Method: equations of motion 

Since the model is noninteracting, its steady-state 
characteristics can be calculated exactly using the 
nonequilibrium Green's function (NEGF) approach 25 . 
This technique has been extensively used in the past 
for studying transport properties in mesoscopic systems 
and molecular junctions 2 ^. We review here the steps 
involved so as to carefully contain the phase factors. 
The derivation presented here follows an equation-of- 
motion approach 2 ^. In this method, an effective quan- 
tum Langevin equation for the subsystem is obtained by 
solving the Heisenberg equations of motion (EOM) for 
the baths variables, then substituting them back into the 
EOM for the subsystem (dots) variables. The indices 
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a,P = 1, 2 identify the two dots. The resulting EOM is 

dap 
~dt 



2 I3 ■ ■ L ■ R 

- -iepap - ir\p - tr)p 



i f drY,^i9i{t-r)C a ,ie^~^ ] a a {r) 

Jtn - i 



10 a,l 



1 \ dT ^2^,r9r (t-r)Ca,r 
J to n.r 



(6) 



The (isolated) reservoirs Green's functions are given by 

g+(t) = -ie-^Oit), g+{t) = -ie-^t). (7) 

The terms vk and rj§ are referred to as noise, induced 
on the subsystem from the left and right reservoirs, re- 
spectively. Their explicit form is 

i 

4 = ^Q,r9r(t~to)a r (t )e-^. (8) 

r 

As an initial condition, we take a factorized state, p(to) = 
PL <8> Pr ® ps(to), with empty dots and the reservoirs 

prepared in a grand canonical state, p v = "ttIp ] ~> 

T v is the temperature of the v — L,R fermi sea and p v 
stands for its chemical potential. The reduced density 
matrix ps denotes the state of the subsystem. Using this 
initial condition, noise correlations satisfy 

i 

(vTitKir)) = E^,re & " (t - T) ^,re i( *?-*"' ) /i i M, 

r 

with the Fermi function f v {u)) = [ e ^-^/ T » + l]^ 1 and 
expectation values evaluated in the Heisenberg represen- 
tation, (A(t)) — Tr[p(to)A(t)]. Steady-state properties 
are reached by taking the limits to — > — oo and t oo. 
We now Fourier transform Eq. ([6]) using the convolution 
theorem with the convention dp(uj) = f x ' oo dtaf3(t)e' lu ' t , 

f]p(uj) — dtr]p(t)e lult . The result, organized in a ma- 
trix form, is 



(9) 



We also define the conjugated-transposed matrix, G~ = 
(G + Y , to be used below. The real part of the self energy 
is a principal value integral, assumed here to vanish. This 
assumption holds when the metals density of states is 
energy independent and the bandwidth is large. We then 
define the hybridization matrix from the relation E + = 
-iT/2, 

T L M = 27re < <**-#> - (12) 



Similar expressions hold for the R side. Using the steady- 
state solution we can write down an expression for 
the reduced density matrix. Back transformed to the 
time domain it takes the form 



(«i«/3) 



Pa,/3 



i r°° r 

-J [(c^G-)^M, 



(G+r K G-' 



a/3 



doj. 



(13) 



The time variable has been suppressed since the result is 
only valid in the steady-state limit. In what follows we 
take £p i and C^,r as rea l constants, independent of the 
level index and the reservoir state, resulting in 



r L _ „ T J{4>p-4>%) V R 



iRd 



(14) 



where the coefficient 7^ , defined through this relation and 
Eq. (|12[) . is taken as a constant, energy independent. 
Using these definitions, the matrix G + takes the form 



G+ = 



^ _|_ i(7-L+7fi) iXLgi0/2 _|_ i2LK e - 



e 

-i<t>/2 



2 

'7-R . 



2 
U! ■ 



■10/2 



£2 



i(~{L+~fn) 
2 



(15) 



and the hybridization matrices are given by 



1L 



1 e^/ 2 

-J0/2 l 



1R 



1 e" 

e i0/2 



1 



We can now calculate, numerically or analytically, the 
behavior of the reduced density matrix under different 
conditions^. Since we are only concerned here with sym- 
metric dot-lead couplings, we take 7^ = 7r = 7/2. Fur- 
thermore, we impose energy degeneracy, ej = £2 = e - 
This choice simplifies the relevant matrices to 



with the Green's function 



The self energies contain the phase factors, 



(10) 



EC^MCcr^--*^- (11) 



G 



u) — e 

22 nn 



¥ feosf 



1 e ^/2 

e -i0/2 2 



1 e -i0/2 
e #/2 I 



(17) 



We present closed analytic expressions for the diagonal 
and off-diagonal elements of the reduced density matrix 
in Sec. III.B. Complementing numerical data for the real- 
time dynamics are included in Sec. IV. This discussion 
is generalized in Sec. V, to include a dephasing probe. 
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B. Observables 

1. Dots occupation 

We expose here two effects that persist away from the 
"symmetric point" , defined as/i^ — e = e — pn'- The dots 
occupations significantly vary with flux, and moreover, 
the degenerate dots acquire different occupations. After 
presenting general expressions away from the symmetric 
point, we consider other relevant cases: the finite-bias 
limit at the symmetric point, the limit of infinite bias 
(which effectively reduces to the symmetric point), and 
the case of 4> = 27rn, n = 0, 1, 2... 

Analytic results are obtained from Eqs. (fT3"|) and (|TT|). 
Organizing these expressions, we find that the occupation 
of dot '1', = (a\ai), is given by two integrals, 



For details, see Appendix A. Since it is a sum of real 
quadratic terms, F± > 0. Inspecting Eq. (f2"Tj) . we note 
that it vanishes in four different cases: (i) at zero bias, 
when pl = (j,r = 0, (ii) at infinite bias, p^ — > oo and 
(Ir • — > — oo, (iii) at the symmetric point when p^ — e = 
e — pn, particularly for e = and ph = — PR, or when (iv) 
<f) = nir, n = 0, 1, 2.. (leading to = F_). Combining 
Eq. (f2"Tj) with the integration of even terms in Eq. JH 
at zero temperature, we resolve the occupations 



J_ 

4-7T 



f L (uj)duj 



(lo — e) 2 + cj 2 — 2lo (lo — e) cos | 



e)] 2 



7 I™ * t \j (^"e) 2 +^o 2 + 2^o(^ 
4t7_oo [(« - e) 2 - w 2 ] 2 + [ 7 (w - e)] 2 

where wc have introduced the short notation 



e)cosf , 

2 , (18) 



ujo — — sm — . 
2 2 



(19) 



- iJ. 



Similarly, the occupation of level '2', p2,2 = {a' 2 a2), is 
given by 



P2,2 



4tt 



f L (uj)duj 



{lo - e) 2 + w 2 + 2w (w - e) cos § 
[( W - e ) 2 - Wo 2 ] 2 + [ 7 ( W 



e)] 2 



"7T / fdu)du] 
47r J-00 [(w 



(cj - e) 2 + cj 2 - 2cj (w — e) cos | 



; ) 2 -c 2 ] 2 + [ 7 ( W - e )] = 



(20) 



In what follows we consider the zero temperature limit. 
The Fermi functions take then the shape of step func- 
tions and the upper limits of the integrals are replaced 
by the corresponding chemical potentials. We now study 
the contribution of the odd term in the integrand. This 
term is responsible for the development of occupation 
difference between the dots, 



_7 

4tt 



- — / diu- 



2wq(w — e) cos § 



8U1 Mn 



[(c- e ) 2 - Wo r + [ 7 (o 



e)] 2 



8tt "~[F- (</>)_ 
where the explicit form of the factors F± is 

■A A 



(21) 



F±(4>) 



2_ 

8 
2 



sin 4 ^ + 2( ML -e) 2 ( Mfl -er 



cos — ± 1 
2 



T C0S 2 Tl 



(Mz 



(MR - e) 2 



(22) 



Pi, 1/2,2 



1 

4tt 
tan" 



27r + tan 



-1 /Mi 



7- 



tan 



Ml 



7+ 



'(^)+--(^) 



smf 

± -In 

8tt 



7- 
F_(0) 



(23) 



The positive sign corresponds to pi.\, the negative sign 
provides p2,2- We have also introduced the short notation 
7± = ^(1± cos|). Equation (j2"31 predicts flux depen- 
dency of electron occupation at degeneracy, using sym- 
metric hybridization constants, once the dots are tuned 
away from the symmetric point. Fig. [2] displays this be- 
havior, and we find that as the dots energies get closer 
to the bias edge, e ~ fiL, the population strongly varies 
with e (panel b). It is also interesting to note that the 
abrupt jump at <f> — 2im (discussed below) disappears 
once the levels reside at or above the bias window, for 
e > ph- This feature results from the strict zero tem- 
perature limit assumed in the analytic calculations. At 
finite T the jump at (f) = 2irn survives even for e > pl- 
However, when the temperature is at the order of the 
hybridization strength, T ~ 7 , the modulation of the 
population with phase is washed out. The following pa- 
rameters are used here and below: flat wide bands, dots 
energies at the order of e = — 0.4, hybridization strength 
7 = 0.05 — 0.5, and a zero temperature, unless otherwise 
specified. The bias voltage is set symmetrically around 
the equilibrium Fermi energy, p^ — —pR, Ap = pl~ PR- 
We now discuss in more details the behavior of the 
occupation in some special cases. First, we consider the 
symmetric point at finite bias and 4> ^ 27rn, n = 0, 1,2.... 
In this case Eq. (|23p precisely reduces to 



Pa,a(^L 



(24) 



This result holds in the infinite bias limit, p^ — > oo and 
pn — oo, irrespective of the (finite) value of e. Next, 
the special case <j> = 2ttti should be separately evaluated. 
At these points we have lo = and Eq. (fT5|) provides 
the simple form at zero temperature 



Pa,a(4> = 27rn) 
7 d.r 

47T 
1 



47T 



oo M ' 7 

_i / PL - e 



1_ 

47T 



dx 



tan 



7 



tan 



Mfl x + T 
-i ( M fl - e 
7 



1 

4' 

(25) 
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FIG. 2: (a) Flux dependency of occupation for dot '1' using 
e = (triangle) e = 0.2 (□) e = 0.3 (o), e = 0.35 (*) and 
e = 0.4 (+). Panel (b) displays results when e is tuned to the 
bias window edge, e ~ /^l, e = 0.29 (□), e = 0.3 (diagonal), 
e = 0.31 (o), and e = 0.31, T = 0.05 (dashed-dotted line). 
In all cases /Ui = — /xr = 0.3, 7 = 0.05, and T = 0, unless 
otherwise stated. 




FIG. 3: (a)-(b) Dots occupations as a function of magnetic 
phase <j) f° r ^-P = 0.6, e = 0.2, T = 0. (c) Occupation 
difference, <5n = pi,i — P2,2. At weak coupling, 7 = 0.05 
(A), the dots occupations are almost identical. When the 
hybridization is made stronger, 7 = 0.5 (o), comparable to 
the levels displacement from the symmetric point, p\ t \ clearly 
deviates from P2,2- At very strong coupling, 7 = 2 (+), the 
occupation difference reduces and asymmetries develop. For 
clarity, results are shown for <f)/-K between (-2,2). 



These points are reflected by abrupt jumps in the 
occupations-flux behavior. Specifically, at the symmetric 
point there is a sharp reduction of occupation number 
from 1/2 [Eq. ([2"3]l] to 1/4 [Eq. ([2"5])]. as observed earlier 
in Ref. 17 . Fig. [5] shows that at strictly zero temperature 
this jump disappears once the dots energies are placed 
at or above the bias edge, e > Thus, the appearance 
of the jump is indicative of the fact that electrons cross 
the junction resonantly. If only tunneling processes con- 
tribute (once the dots' energies are placed above the bias 
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FIG. 4: Occupation difference as a function of bias voltage 
Ap, for different magnetic flux values, (j> = n/2 (full line), 
4> = 7r/4 (dashed line), <f) — n/8 (dashed-dotted line). Other 
parameters are e = 0.2 and 7 = 0.05, T = 0. The inset 
presents data for backward and forward biases; the main plot 
zooms on the positive bias regime. 
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FIG. 5: Occupation difference as a function of magnetic flux 
for different bias values, Afi — 0.1 (full line), A/x = 0.2 
(dashed line), A^i — 0.3 (dashed-dotted line) and Afi — 0.4 
(dotted line). Other parameters are e = 0.2, 7 = 0.05 and 
T = 0. 



window and the temperature is very low) , the population 
continuously vary with flux. 

The total electronic occupation of the dots, at steady- 
state, generalizes the standard symmetric case attained 
in RefA 



Pl,i + P2.2 



J_ 

2tt 



duj 
00 [(u) 



[{uj-ef+ul] [f L {u) + f R {u)] 



(26) 



We now highlight one of the main results of the paper, the 
onset of occupation difference in this degenerate (ej = e 2 ) 
and spatially symmetric (7^ = 7^) setup. Using Eq. 
(|23l), we find that 



5n 



- P2,2 



-in 



In 



(27) 



As we mentioned above, this quantity is nonzero when 
the following (sufficient) conditions are simultaneously 
satisfied: (i) the bias voltage is finite, neither zero nor 
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is only determined by the offset of e from the center of 
the bias window. In Fig. [5] we further show the flux 
dependency of Sn, which is particularly significant when 
Ap = 2e. 

The effect of finite temperature on the occupation-flux 
dependence, and on the development of occupation dif- 
ference, is displayed in Fig. [BJ We find that the effects 
largely survive at finite T, as long as T < 7. These re- 
sults were calculated numerically, based on Eqs. (IT8l) and 



FIG. 6: Finite temperature effect. Main plot: Occupation 
difference as a function of bias voltage for (j> = n/4. Inset: 
occupation of dot '1' as a function of magnetic phase for A/j, = 
0.4. In both panels T = (dotted line), T = 0.01 (full line), 
T = 0.05 (dashed-dotted line) and T — 0.1 (dashed line). 
Dots parameters are e = 0.2 and 7 = 0.05. 

infinite, (ii) the dots are positioned away from the sym- 
metric point, e 7^ +[Xft)/2, and (ii) the phase <f> is not 
a multiple of tt, <f> 7^ rm, n = 0, 1, 2.... To rephrase this 
observation, the occupation difference can be controlled 
by manipulating the subsystem-metal hybridization en- 
ergy 7, by changing the bias voltage, by applying a gate 
voltage for tuning the dots energies, and by modulating 
the phase 4> through the magnetic flux. The role of these 
control knobs are illustrated in Figures [3J 0] and [5] 

In Fig. [3] we display the levels occupation in the reso- 
nant regime, p,R < e < /i^ while varying 7. At weak cou- 
pling Sn is insignificant. However, the occupation differ- 
ence becomes large when co-tunneling effects contribute. 
More notably, Fig. @] illustrates the strong controllability 
of Sn with applied voltage. We find that the occupation 
difference is maximized at the edge of the resonant trans- 
mission window, when \il — e = (or equivalently, when 
Afi = 2e). The magnetic phase affects the width and 
height of the peak, but not the absolute position which 



2. Coherence 

It was recently argued that the decoherence behav- 
ior in our generic setup, including two noninteracting 
(uncoupled) quantum dots interferometer, can be sup- 
pressed when the device geometry is made asymmetric 
and nondegenerate, using e\ 7^ £2 and 7^ 7^ 7a 19 . The 
requirement for asymmetry in this work arises from the 
observation of the "phase localization" effect, which hin- 
ders phase manipulation in the system at the symmet- 
ric point. The term "phase localization" refers to the 
fact that if we define pi^(t) — |pi,2(i)|e , *'^' ) the relative 
phase ip localizes to the values —it/ 2 or 7r/2 in the long 
time limit when (f> 7^ 27m, n is an integer—. Based on nu- 
merical simulations, we have pointed out in Ref. 20 that 
phase localization occurs only at the symmetric point, 
while at other values of e the real part of pi t 2 is finite 
and nonzero in the asymptotic limit for any phase be- 
sides 2irn—. This observation is established here analyti- 
cally in the steady-state limit, implying that decoherence 
could be suppressed in degenerate-symmetric systems by 
gating the dots, shifting their energies relative to the bias 
window. 

We derive a closed expression for the off-diagonal sys- 
tem element pi.2 = (a|&2) by studying Eq. (|13p . 



II"" I cos 

P12 = / fL{u)duj 



I [(to - e) 2 - u 2 ] +zsin | [( LU -e) 2 +uj l 



4tt 



f R (uj)duj 



[( W - e ) 2 -c 2 ] 2 + [7(^-6)] 2 J 
cos I [{w - e) 2 - oj 2 ] - i sin f [{u - e) 2 + u 2 ] 
[(^- e ) 2 -w 2 ] 2 + [7(^- £ )] 2 

I 



(28) 



At finite bias and zero temperature direct integration 
provides the real (JR.) and imaginary (3) parts of p\^. 
(cf) ^ 2im) 



and 



3/3i,2 = -j- sin(</>/2) 



tan 



-1 ( VL ~ e 

7+ 



tan 



-1 (Vl - e 



7- 



7+ 



7- 



tarT 



(30) 



tan 



i(^)_ tan -i(^p)] ; (29) 



As before, we define j± = ?(1 ± cos t). We now readily 
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FIG. 7: Real and imaginary parts of the coherence as a func- 
tion of the bias voltage. <f> — n (dashed line), <j> — ir/2 (full 
line), <j> — 7r/4 (dashed-dotted line). Other parameters are 
e = 0.2, 7 = 0.05, and T = 0. The oval shape marks the 
region of phase localization at positive bias. 
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FIG. 8: Transmission coefficient as a function of energy for 
e = 0.2, 7=0.05, = tt (dotted line), <f> = tt/2 (dashed line) 
and (j) — (full line). The inset presents the charge current 
for <j> = 7r/2 and e = 0.2 (dashed line), e = (dotted line). 



obtained by defining the number operator Nl = ajai, 
providing the current Jl->r = —^f- = —i[H,Nz]. This 
yields 



confirm that at the symmetric point the real part van- 
ishes and "phase localization" takes place 17 . In particu- 
lar, in the infinite bias limit we find 3pi.2 = ^sin^, in 
agreement with previous studies^. We also include the 
behavior at the special points <f> = 2irn. Eq. (|28|) reduces 
then to a simple Lorentzian form, at zero temperature, 



Pl.2\ 



47T 



dx 



1_ 

4tt 



dx 



1 

47T 



tan 



PL - e 



x 2 + 7 2 

-l / P<; - < 



tan 



7 



(31) 



The sign reverses for <f> — ±2tt. We note that the imag- 
inary part of the coherence identically vanishes at zero 
phase while the real part is finite, approaching the value 
1/4 at the symmetric point. 

Numerical results in the steady-state limit are dis- 
played in Fig. [JJ We find that both the real and imagi- 
nary parts of /0x,2 demonstrate significant features when 
the dots' levels cross the bias window, at Afi = 2e. The 
value of the real part abruptly changes sign, the imagi- 
nary part develops a step. At large bias Sft/Oi^ diminishes 
while 9/9i,2 is finite, indicating on the development of 
the phase localization behavior. It can be shown that 
the double-step structure of 3/01,2 (as a function of A/i) 
disappears when the dots energies are set at the symmet- 
ric point. 



3. Current 

It is of interest to complement our study of subsystem 
(dots) properties and examine the transmission coeffi- 
cient and the overall charge current in the system. The 
electric current, flowing from the L metal to the R end, is 



2,a=l,2 



{a\a a ) 



(32) 



Expectation values are calculated in the steady-state 
limit. Using the EOM formalism as explained in Sec. 
Ill, we get the standard result 



Jl- 



1 

2^ 



dLuT LR (Lu)[f L (u;) - f R (iu)}. (33) 



The transmission coefficient is defined as Tlr = 
Ti^r^G+r^G - ), where the trace is performed over the 
states of the subsystem (dots). In the present model, at 
zero temperature, we obtain 



Jl- 



■R 



1 



) 2 cos 2 % 



-e) 2 



2l 2 



<t> 

COS 77 

2tt 



7 



I V 7+ 



_ 7 _ {ta „-.(^ 



J 2 (uj — e) 2 

--(^)} 



\ 7_ 



■)}], (34) 



which agrees with known results^. Using the NEGF for- 
malism, we could similarly investigate the shot noise in 
the double-dot AB interferometer—. 

The transmission function is plotted in Fig. [5] dis- 
playing destructive interference pattern for (f> — 7r and a 
constructive behavior for = 0. For </> ^ mr the trans- 
mission nullifies exactly at the position of the resonant 
levels. The inset presents the current-voltage character- 
istics for (p — 7r/2 away from the symmetric point (dashed 
line), and at the symmetric point (dotted line). We note 
that the double-step structure disappears at the latter 
case. It can be shown that the double step structure of 
(see Fig. [JJ similarly diminishes at the symmetric 

point. 
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IV. TRANSIENT BEHAVIOR 

It is of interest to investigate the development of the 
phase dependency of the occupancy, and the occupancy 
difference Sn, before steady-state sets. Similarly, the dy- 
namics of coherences is nontrivial even without electron- 
electron interaction effects 2 ^. We complement the NEGF 
steady-state expressions of Sec. Ill with numerical calcu- 
lations of the transient behavior using an exact numerical 
tool that is based on the fermionic trace formula^ 



Tr [ 



e M lg M 2 



2 ...e Mp ] =det[l 



(35) 



Here m p is a single-particle operator corresponding to a 



quadratic operator M p = X^j( to p 



! (oj) are 

fermionic creation (annihilation) operators. The trace is 
performed over all electronic degrees of freedom. Our 
objective is the dynamics of a quadratic operator B = 
1,2, 



a\a k , j, k 



(B(t)) = Tr [p(t o y m Be- im ] 

= \mi^^Tr[p L p R p s e im e XB e- lHt } .(36) 

We introduce the A parameter, taken to vanish at the 
end of the calculation. The initial condition is factorized, 
p(ta) — ps ® PL ® PR, and these density operators follow 
an exponential form, e M , with M a quadratic operator. 
The application of the trace formula leads to 



(e XB ^) = det [\I L - f L ] ® [I R - f R ] ® [I s - f s ] 

7z®/n®/s}- (37) 



with b and h as the single-body matrices of the B and 
H operators, respectively. The matrices l v and 7s are 
the identity matrices for the v = L,R space and for the 
subsystem (dots). The functions Jl and fn are the band 
electrons occupancy /„(e) — [e^ e_AI ' y ) + 1] _1 . Here they 
are written in matrix form and in the energy represen- 
tation, fs represents the initial occupation for the dots, 
assumed empty, again written in a matrix form. When 
working with finite-size reservoirs, Eq. (|37[) can be read- 
ily simulated numerically-exactly. 

Fig. IH1 displays the evolution of the occupation differ- 
ence, presented as a function of A/i. In this simulation 
we used finite bands with a sharp cutoff, D = ±1. At 
short time Sn shows weak sensitivity to the actual bias. 
Only after a certain time, jt ~ 2, the peak around the 
edge at Ap = 2e clearly develops. Note that since the 
band is not very broad, edge effects are reflected at large 
biases as nonzero occupation difference, in contrast to 
the broad-bandwidth long-time behavior of Fig. @] 

The transient behavior of the coherences, Sftpi^ and 
Sypi,2> is included in Fig. QJJ the corresponding steady- 
state value are presented in Fig. [7] We can follow the 
temporal features of the phase localization effect, i.e., the 
disappearance of the real part of the coherence at the 
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FIG. 9: Time evolution of the occupation difference, 7 = 0.05, 
e = 0.2, T = 5 x 10 _s , </> = ty/2. 
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FIG. 10: Time evolution of the real (a) and imaginary (b) 
parts of the coherence 7 = 0.05, e = 0.2, T — 5 x 10~ 3 , 
(j> = tt/2. 



symmetric point or at large bias, when <\> 7^ 27m. Using 
<f> = tt/2 we note that while at short to intermediate 
time (jt < 2) significant coherence builds up, the real 
part of the coherence eventually survives only at small 
biases. Regarding timescales, we find that while 3/31,2 
reaches the steady-state values at short time, jt ~ 2, 
5£/?i,2 approaches its stationary limit only at longer times, 
for jt ~ 10. Similar results were obtained in Ref^S. 



V. DEPHASING PROBE: STEADY STATE 
CHARACTERISTICS 

We have discussed so far pure coherent evolution ef- 
fects in double-dot AB interferometers. It is important 
to examine at this point the role of elastic dephasing ef- 
fects on this evolution as was done experimentall y 3 ^ 34 
and theoretically 3 ^— in related systems. Here, we are 
essentially focused on the effect of dephasing on the mod- 
ulation of occupation with magnetic phase. Recently, 
interference and decoherence processes were studied not 
only in quantum dot structure o 23 ' 24 , but in molecular 
loops as welt^ 1 ^^. 

Phase-breaking processes arise due to the interaction 
of electrons with other degrees of freedom, e.g., with elec- 
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trons, phonons and defects. We generalize here the dis- 
cussion of Sec. Ill, and incorporate dephasing processes 
into our system phenomenologically, by using the well 
established method of Biittiker dephasing probed. In 
this technique, elastic dephasing processes on the dots 
are emulated by including a third terminal, P, enforc- 
ing the requirement that the charge current towards the 
probe terminal, at a given electron energy, should vanish. 
Thus, electrons travel to the probe and return to the sys- 
tem with a different phase, while both electron number 
and electron energy are conserved. This condition sets 
an electron distribution within the probe. As we show 
below, away from the symmetric point this distribution 
effectively depends on the magnetic phase. Other phe- 
nomenological tools to incorporate dephasing processes 
in mesoscopic devices are based on the introduction of 
random-phase fluctuations into the scattering matrix-^, 
or on the inclusion of damping terms into the off-diagonal 
elements of the density matrix within quantum master 
equations (Lindblad or Redfield) formalisms^. 

Using Biittiker probe method, we augment the Hamil- 
tonian ((T|) with a probe, adding to the system a nonin- 
teracting electron reservoir P, 



Hn = H 



P&P 



\ p a\a p + h.< 



(38) 



oeP 



The parameter A denotes the coupling strength of dot '1' 
to the P terminal, taken as a real number. Note that we 
only allow here for local dephasing on dot T'. One could 
similarly consider models where both dots are suscepti- 
ble to dephasing effects, possibly from different sources. 
Following the equations-of-motion approach as detailed 
in Sec. Ill, we arrive at the steady-state expression for 
the reduced density matrix 



,t, 



1 f°° 

Zj " . r D D J — CO 



v=L.R.P 



The probe hybridization matrix is given by 



r p = 7 p 



1 o 





(39) 



(40) 



and the dot's Green's function is written by generalizing 
the matrix (fl"5)) . to include the probe self energy, 



G+ = 



IT 9 9. '9 



-i<t>/2 



+ 



2 



i>/2 



-i<j>/2 



€ 2 H 2 



This matrix is written here in a general form, to allow 
one to distinguish between the two dots and the differ- 
ent dots-metals hybridization terms. The dot-probe hy- 
bridization is defined as = 2tt^2 p |A p | 2 <5(u; — wo), in 
analogy with Eq. (|12[) . In our calculations below we as- 
sume energy degenerate dots and symmetric couplings, 
e = ei = e 2 , jl = Jr = 7/ 2 - 

We now derive the probe distribution by demanding 
that the energy resolved charge current to the P terminal 
vanishes. The total current to P is given by the sum of 
the currents from the L and R terminals, generalizing 
Eq. (33), 



Jl->P + JR-yP 
1 [°° 

— \ Tlp{u) [/l(w) - /pH] duj 



J P 



2tt 



Tnp(u)[fR(u)-f P (u)]du (41) 



with the transmission coefficient T u v{u) = 
Trp^G+r^G - ]. By requiring the integrand to vanish, 
we arrive at the probe distribution 



/pM 



7lp(^)/l(^) +TRp{uj)fR.{uj) 
Tlp(u) + Trp(w) 



(42) 



Direct evaluation of these transmission coefficients pro- 
vide the electron distribution in the probe, 



j(u — e) sin % cos % r , , , „ , , 
" 2 [/lH-ZaH]. (43) 



2[( W -c)a+ w g] 



As before, u>o = ^sin^. This expression indicates that 
the magnetic flux plays a role in setting the distribution 
within the probe, (such that it only dephases the system 
and does not deplete electrons or allow energy reorga- 
nization). This dependency disappears when the dots 
energies are set at the symmetric point, since the contri- 
bution of the second term in Eq. (|43| diminishes in the 
integrals of Eq. (|44|) . from symmetry considerations. We 
now write integral expressions for the dots occupations 
using Eq. ([3"9"1) . 
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2tt, 
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(w - e) 2 + cj 2 + 2w (w - e) cos ^ 



(cj — e) 2 + w 2 + 2w (w — e) cos — + w 7p sin ^ + 



f L (oj) 



7p 



(u) — e) + w - 2w (<^ — e) cos — + w 7p sin — H — — 



/flM 



7 IP 2 <P 

-t; — cos 77 
8tt 2 



d,LU 



/p(w), 



(44) 



with 



A(w) = (o> - e) 



2 77P , • / 7p\, s 



In the absence of dephasing these expressions reduce to 
Eqs. (fT8|) and ([20]) . In the opposite limit, at very large 
dephasing, jp 7, jp > A/i, we note that p2,2 is dom- 
inated by 7p7 terms that are flux independent, while 
is dominated by its last term, oc "fpfp, which is 
flux dependent away from the symmetric point, result- 
ing in pi i cx sin(0). Thus, quite counter-intuitively, we 
find that the level that is directly susceptible to local de- 
phasing demonstrates flux dependency of occupation at 
strong dephasing, while the level that indirectly suffers 
dephasing effects more feasibly looses its coherent oscil- 
lations. 

Using numerical integration, dots occupations and 
their oscillation with phase are presented in Fig. [TTJ Wc 
observe the following trends upon increasing dephasing 
strength: At the symmetric point, (a)-(b), the abrupt 
jump at zero magnetic phase immediately disappears 
with the application of finite dephasing. When the dot 
energies are placed away from the symmetric point, yet 
they buried within the bias window, (c)-(d), the abrupt 
jump at zero magnetic phase again disappears, though 
the oscillations of occupation with phase prevail till large 
dephasing, jp ~ A/i. More significantly, when the dots 
energies are tuned at the edge of the bias window, (e)-(f), 
we find that dot T' (which is directly dephased) devel- 
ops new type of oscillation with phase. Only at very large 
dephasing, "fp 3> A/Lt, these oscillations are overly sup- 
pressed. Thus, away from the symmetric point not only 
features of coherent dynamics survive even at significant 
dephasing strength, new type of coherent oscillations may 
develop as a result of the application of clastic scatter- 
ing effects on the dots. It is interesting to reproduce this 
behavior while modeling elastic dephasing effects using 
other techniquesj^ 2 .^— . 
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FIG. 11: The role of dephasing on the dots occupations - 
magnetic phase dependency, (a)-(b) e = 0, (c)-(d) e = 0.2, 
(e)-(f) e = 0.3, where (dotted line), jp = 0.01 (dashed 

line), 7p = 0.05 (dashed-dotted line) and fp — 0.5 (full line). 
Other parameters are 7 = 0.05, A/j, — 0.6, T = 0. 



VI. CONCLUSIONS 

In this paper, we have addressed the issue of magnetic 
field control on electronic occupation and coherence in 
double-dot AB interferometers. The system under inves- 
tigation included energy degenerate dots with symmetric 
dot-metals hybridization strengths. However, by voltage 
gating the dots energies away from the so-called sym- 
metric point at which e = (fip + /j,r)/2 we have resolved 
four nontrivial effects that can allow for significant con- 
trollability over dots occupations and their coherence: 
(i) Dots occupations may significantly vary with mag- 
netic flux, particularly when the dots levels reside close 
to the bias edge, (ii) The dots acquire different occupa- 
tions, though they are energy degenerate. This behav- 
ior is maximized at the bias edge e ~ /xj,. It survives 
at finite temperature, as long as T < 7. (iii) Regard- 
ing the dots coherence, we have proven that the effect 
of "phase localization"— does not take place away from 
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the symmetric point, allowing for decoherence control in 
the system. Furthermore, (iv) we have found that away 
from the symmetric point the system can withstand de- 
phasing processes, maintaining its coherent evolution and 
even develop new type of oscillations under intermediate 
dephasing strengths (~fp ~ A/x, and 7p ^> 7). 

Our minimal model could be applied to describe mag- 
netic field control in mesoscopic conducting loops and in 
molecular ring structures. In the latter case it was partic- 
ularly noted that degeneracy is crucial for allowing con- 
trollability within realistic magnetic fields^. Our study 
has been limited to the noninteracting electron model, 
excluding both electron-electron interaction effects and 
other explicit sources for dephasing and inelastic scat- 
tering processes. It is of interest to explore the role of 
interactions on the effects revealed in this paper, as we 
expect it to lift the energy degeneracy in the system, 
further intensifying the effects discussed here. This be- 
havior can be immediately observed at the mean-field 
level. The Hartree term corrects the dots energies, e.g. 
(\ (\ + Up2,2r^- Thus, away from the symmetric point, 
(flux generated) unequal dots occupations translate to ef- 
fective unequal energetics for the two (identical) dots. 
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where d = (fin — e), c — (fiL — e), b — 7 and a = \ sin |, 
leading to by/Aa 2 — b 2 — ±ij 2 cos ^ and 



2a 2 - b 2 - 2d 2 = 7 2 



2a 2 - b 2 - 2c 2 = 7 2 



- sin 1 



1.2* , 

-sin — — 1 



- 2(Mfl - ef 



2(ML ~ e) 2 



2 2 

We now reorganize Eq. (|A2I) using the relations 
tan -1 a; + tan -1 y = tan -1 ( \-x y J an< ^ tan -1 z 
\ [ln(l - iz) - ln(l + iz% to find 
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"(Mi - e) 2 W sin 2 77 - (fi R - e) 2 - 7 2 ( 1 ± cos ~ 
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(A4) 



We can also reorganize these factors as a sum of real 
quadratic terms, 



Appendix A: Derivation of Equation (|21|) 

In this appendix our aim is to evaluate the following 
integral analytically 
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, m [(u-e) 2 -oj 2 r + [^-e)}- 
We achieve this by using the following definite integral 
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Attaching the missing prefactors, / = j^2wocos^/o, we 
obtain Eq. (J2lj) 



8tt 



In 



F-{<t>) 



(A6) 



1 G. Hackenbroich, Phys. Rep. 343, 463 (2001), and refer- 
ences therein. 

2 Y. Imry, Introduction to Mesoscopic Physics, 2nd ed. (Ox- 
ford University Press, Oxford, 2002). 

3 A. Yacoby, M. Heiblum, D. Mahalu, and H. Shtrikman, 
Phys. Rev. Lett. 74, 4047 (1995). 

4 R. Schuster, E. Buks, M. Heiblum, D. Mahalu, V. Uman- 
sky, and H. Shtrikman, ature 385, 417 (1997). 



5 J. Konig and Y. Gefen, Phys. Rev. Lett. 86, 3855 (2001); 
Phys. Rev. B 65, 045316 (2002). 

6 F. Li, X.-Q. Li, W.-M. Zhang, and S. A. Gurvitz, Euro. 
Phys. Lett. 88, 37001 (2009). 

7 F. Li, H.-J. Jiao, J.-Y. Luo, X.-Q. Li, S. A. Gurvitz, Phys- 
ica E 41, 1707, (2009). 

8 Y. Tokura, H. Nakano, and T. Kubo, New J. Phys. 9, 113 
(2007). 



12 



9 Y.-S. Liu, H. Chen, and X.-F. Yang, J. Phys.: Condens. 
Matt. 19, 246201 (2007). 

10 D. Boese, W. Hofstetter, and H. Schoeller, Phys. Rev. B 
66, 125315 (2002). 

11 V. Kashcheyevs. A. Schiller, A. Aharony, and O. Entin- 
Wohlman, Phys. Rev. B. 75, 115313 (2007). 

12 A. Aharony and O. Entin-Wohlman, Phys. Rev. B 72, 
073311 (2005). 

13 O. Hod, E. Rabani, and R. Baer, Ace. Chem. Res. 39, 109 
(2006). 

14 O. Entin-Wohlman and A. Aharony, Phys. Rev. B 85, 
085401 (2012). 

15 D. Rai, O. Hod, and A. Nitzan, J. Phys. Chem. Lett. 2, 
2118 (2011), Phys. Rev. B 85, 155440 (2012). 

16 D. Rai and M. Galperin, Phys. Rev. B 86, 045420 (2012). 

17 M. W.-Y. Tu, W.-M. Zhang, and J. Jin, Phys. Rev. B 83, 
115318 (2011). 

18 M. W.-Y. Tu, W.-M. Zhang, J. Jin, O. Entin-Wohlman, 
and A. Aharony, Phys. Rev. B 86, 11543 (2012). 

19 M. W.-Y. Tu, W.-M. Zhang, and F. Nori. larXiv:1204.5931l 

20 S. Bedkihal and D. Segal, Phys. Rev. B 85, 155324 (2012). 

21 B. Kubala and J. Konig, Phys. Rev. B 65, 245301 (2002). 

22 M. Biittiker, Phys. Rev. B 32, 1846 (1985); Phys. Rev. B 
33, 3020 (1986). 

23 S. Pilgram, P. Samuelsson, H. Forster, and M. Biittiker, 
Phys. Rev. Lett. 97, 066801 (2006). 

24 H. Forster, P. Samuelsson, S. Pilgram, and M. Biittiker, 
Phys. Rev. B 75, 035340 (2007). 

25 Y. Meir and N. Wingreen, Phys. Rev. B 68, 2512 (1992). 

26 J. Franson, Non- equilibrium nano-physics, a many body ap- 
proach, Lecture Notes in Physics 809, Springer (2010). 

27 A. Dhar and D. Sen, Phys. Rev. B 73, 085119 (2006). 



28 M. Biittiker, Phys. Rev. B 46, 12485 (1992). 

29 A. Batra, G. Kladnik, H. Vzquez, J.S. Meisner, L. Flore- 
ano, C. Nuckolls, D. Cvetko, A. Morgante, L. Venkatara- 
man, Nature Comm. 3, 1086, (2012). 

30 I. Klich, in "Quantum Noise in Mesoscopic Systems", 
edited by Yu. V. Nazarov and Ya. M. Blanter (Kluwer, 
2003). 

31 M. G. Pala and G. Iannaccone, Phys. Rev. B 69, 235304 
(2004). 

32 H.-P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press, Oxford, UK, 
2002). 

33 D. Rohrlich, O. Zarchin, M. Heiblum, D. Mahalu, and V. 
Umansky, Phys. Rev. Lett. 98, 096803 (2007). 

34 D.-I. Chang, G. L. Khym, K. Kang, Y. Chung, H.-J. Lee, 
M. Seo, M. Heiblum, D. Mahalu, and V. Umansky, Nature 
Physics 4, 205 (2008). 

35 A. Ueda and M. Eto, Phys. Rev. B 73, 235353 (2006); New 
J. of Phys. 9, 119 (2007); Physica E 40, 1602 (2008). 

36 Z. Zhu, A. Aharony, O. Entin-Wohlman, and P. C. E. 
Stamp, Phys. Rev. A 81, 062127 (2010). 

37 T. Kubo, Y. Tokura and S. Tarucha, J. Phys. A: Math. 
Theor. 43 354020 (2010). 

38 R. Haxtle, M. Butzin, O. Rubio-Pons, and M. Thoss, Phys. 
Rev. Lett. 107, 046802 (2011). 

39 S. Ballmann, R. Hartle, P. B. Coto, M. Elbing, M. Mayor, 
M. R. Bryce, M. Thoss, and H. B. Weber, Phys. Rev. Lett. 
109, 056801 (2012). 

40 M. Sindel, A. Silva, Y. Oreg and J. von Delft, Phys. Rev. 
B 72, 125316 (2005). 



